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We present updated results on time-dependent CP asymmetries in fully reconstructed 
B _ > .D?*) ± 7r = F alm B°-^D ± p T decays in approximately 232 million T(4:S)-*BB events collected 
with the BABAR detector at the PEP-II asymmetric-energy B factory at SLAC. From a time- 
dependent maximum likelihood fit we obtain for the parameters related to the CP violation angle 
2/3+ 7 : 

a Dw = -0.010 ± 0.023 ±0.007 , cg£ = -0.033 ± 0.042 ±0.012 , 
a"' 1 = -0.040 ± 0.023 ±0.010 , cg*^ = 0.049 ± 0.042 ±0.015 , 
a Dp = -0.024 ± 0.031 ± 0.009 , = -0.098 ± 0.055 ± 0.018 , 

where the first error is statistical and the second is systematic. Using other measurements and 
theoretical assumptions, we interpret the results in terms of the angles of the Cabibbo-Kobayashi- 
Maskawa unitarity triangle and find | sin(2,3+7)j >0.64 (0.40) at 68% (90%) confidence level. 

PACS numbers: 12.15.Hh, 11.30.Er, 13.25.Hw 



In the Standard Model, CP violation in the weak in- 
teractions between quarks manifests itself as a non-zero 
area of the Cabibbo-Kobayashi-Maskawa (CKM) unitar- 
ity triangle [ij. While the measurement of sin 2/3 is now 
quite precise 0, > the constraints on the other two an- 
gles of the unitarity triangle, a and 7, are still limited by 
statistical and theoretical uncertainties. 

This paper presents updates for the measurements of 
CP asymmetries in B°^D^ ± n zf: decays [4J, as reported 
in Ref. [5j, with a larger data sample (x2.6), and in addi- 
tion includes the measurement of the CP asymmetry in 
the decay mode B°~^D ± p T . We denote these decays as 

is a charged pion or p meson. 

The time evolution of B°~^D^ ± h T decays is sen- 
sitive to 7 because the CKM-favored decay amplitude 
B a ^D^ + h~ , which is proportional to the CKM ma- 
trix elements V cb V* d , and the doubly-CKM-suppressed 
decay amplitude B°— >£>(*) + /i~ , which is proportional to 
V cd V* b , interfere due to B° - B° mixing. The relative 
weak phase between these two amplitudes is 7. With 
B° — B° mixing, the total weak phase difference between 
the interfering amplitudes is 2/3+7. 

Neglecting the very small decay width difference be- 
tween the two B° mass eigenstates |6j, the proper-time 
distribution of the B°^D ( -^ ± h T decays is given by 



Ai) 



3 -|At|/ T 

4r 



x [1 + S c sin(Am d At) (1) 



+ rj C cos(Am d Ai)] , 

where r is the B° lifetime, Am^ is the B Q —B° mixing 
frequency, and At = t rec — t tag is the time difference 
between the B°-^D^ ± h T decay (B rec ) and the decay of 
the other B (Stag) from the T(4S) -> B°B° decay. In this 
equation the upper (lower) sign refers to the flavor of B tag 
as B° (B°), while r, = +1 (-1) and C = ± (-) for the 



final state D^-fi + (D^ + h~). The sine term is due to 
interference between direct decay and decay after B°— B° 
mixing. The cosine term arises from interference between 
decay amplitudes with different weak and strong phases 
(direct CP violation) or from CP violation in mixing. 
The S and C asymmetry parameters can be expressed as 



S± = - 



2Im(A±) 
1 + |A±| 2 



and C = 



1 - r A 



where r = |A_i 



1/IAJ and 



g AjB^D^h*) 
p A{B°-^D(*^h ± ) 



r ±l e -i(2/3+7T<5\ 



(2) 



(3) 



Here - is a function of the elements of the mixing Hamil- 
tonian 6] , and <5 is the relative strong phase between the 
two contributing amplitudes. In the Standard Model, 
CP violation in mixing is negligible and thus |^| = 1. In 
these equations, the parameters r and 6 depend on the 
choice of the final state. They will be indicated as r D7T , 



8°* for the B°^D ± tt^ 
and r D * v , 5 D ** for B°- 



mode 



„d p s d p for B o 



HI 



>D±pT, 



Interpreting the S parameters in terms of the an- 
gles of the unitarity triangle requires knowledge of the 
corresponding r parameters. The values of r are ex- 
pected to be small (^0.02) and therefore cannot be ex- 
tracted from the measurement of C. They can be esti- 
mated, assuming SU (3) symmetry and neglecting con- 
tributions from Il^-exchange diagrams, from the ratios of 
branching fractions B(B ^D ( * )+ n-)/B(B°^D^-ir+) 
and B(B°-+Df p-)/B{B a ^D- p+) UuME^ 

This measurement is based on 232 million TJAS) — > 
BB decays, collected with the BABAR detector |l2j at the 
PEP-II asymmetric-energy B factory at SLAC. We use a 
Monte Carlo simulation of the BABAR detector based on 



GEANT4 [13j to validate the analysis procedure and to 
estimate some of the backgrounds. 

The event selection criteria are unchanged from our 
previous publication [jj, except for the application of 
a kaon veto on the pion candidate in the decay modes 
£)(*)- 7r + to suppress B°—*D(*)~K~^ background events, 
and for the addition of the decay mode B°—>D~p + . The 
D*~ is reconstructed through its decay to D°tt~ , where 
the D° decays into K + ir~, K + ir~ir°, K + ir~ir + ir~ , or 
K®ir + n~. The D~ is reconstructed through its decay 
into K + ir~TT~ or K®ir~ 



in the final state tt +i 
the 7r + 7T° invariant mass 



The p + decay is reconstructed 
For the CP analysis we require 
m^K ) to be in the window 
620 < m W7V o < 920 MeV/c 2 . Exploiting the polarization 
of the p meson from the decay B°^D~ p + we require the 
cosine of the p + helicity angle 6*h e i, defined as the angle 
between the charged pion and the D~ momentum in the 
p + rest frame, to satisfy | cos^hcil >0.4. 

The beam-energy substituted mass, toes = 
yj s/4 — p* B 2 , and the difference between the B 
candidate's measured energy and the beam energy, 
AE = Eg - (s/s/2), are used to identify the final 
sample, where E* B (p* B ) is the energy (momentum) of 
the B candidate in the nominal e + e _ center-of-mass 
frame, and t/s is the total center-of-mass energy. The 
AE signal region is defined as \AE\ < 3a, where the 
resolution a is mode-dependent and approximately 
20 MeV, as determined from data. Figure shows the 
toes distribution for candidates with toes > 5.2 GeV/c 2 
in the AE signal region. These candidates satisfy the 
tagging and vertexing requirements, which are described 
later. Each distribution is fit to the sum of a thresh- 
old function [l4|, which accounts for the background 
from random combinations of tracks (combinatorial 
background), and a Gaussian distribution with a fitted 
width of about 2.5 MeV/c 2 , which describes the signal 
and the backgrounds that peak in the toes signal 
region (peaking background). Signal yields and sample 
purities are determined in the toes signal region, with 
toes > 5.27 GeV/c 2 , and are summarized in Tabled 
Backgrounds from B° and B + decays that peak in the 
TOgg signal region are estimated using Monte Carlo 
events and are mostly due to charmed final states. They 
are also reported in Table [I] 



TABLE I: Signal yields, sample purities P, and fractions 
of peaking backgrounds, f pea k, for the selected samples for 
events that satisfy the tagging and vertexing requirements 
described in the text. 



Decay 
mode 


Yield 


P(%) 


fpea 

B° 


c(%) 

B ± 


B~*D ± p" F 


15038 ± 132 
14002 ± 123 
8736 ± 101 


87.0 ±0.3 
93.2 ± 0.2 
81.7 ±0.4 


1.6 ±0.1 
1.0 ±0.1 
1.3 ±0.2 


1.2 ±0.1 
1.1 ±0.1 
1.5 ±0.2 
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B -> Dlt 

Purity: 87.0 ±0.3% 
Yield: 15038 ±132 




B -> D It 

Purity: 93.2 ±0.2% 
Yield: 14002 ±123 
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Purity: 81 .7 ±0.4% 
Yield: 8736 ±101 
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FIG. 1: tties distributions in the signal region for, from top to 
bottom, the B^D ± -r t , B-tD****, and B->D ± p* sample 
for the events that satisfy the tagging and vertexing require- 
ments described in the text, fit with the function described in 
the text. The dashed lines indicate the sum of the combina- 
torial and peaking background contributions. 



For the B°^D ± p T mode we consider additional 
sources of background with the same final state _D ± 7r =F 7r°, 
where the tt^tt system is not produced through the p^ 
resonance. Interfering sources of background can intro- 
duce a dependence of the A^ p parameters of Eq. [3] on 
m W7r o. The dependency has been studied using the dis- 
tribution Of TOjr^o . 

The possible background contributions have been eval- 
uated with a sample of 130273 B°^D~ir + ir° candidates, 
on which the requirements on the p helicity and on m nv o 
have been removed. Three interfering components are 
considered: B°^D-p+ (the signal), B°^D- p'+ (1450) 
with a pole mass of (1465 ± 25) MeV/c 2 and a width of 
(400 ± 60) MeV/c 2 @ for the p'+, both described with P- 
wave relativistic Breit-Wigner functions [l5L fl6| , and a 
non-resonant component, B°— >D~ (7r + 7r°)„ r . Contribu- 
tions from the decay modes B°^D*~ir + (D*~ -^D~ir°) 
and B°^D**°ir° (D**°^D-tt+) are negligible due to the 
kinematic constraints imposed on the p daughter parti- 
cles. We perform a fit to the binned to^o distribution to 
extract the amplitudes of the three components, where 
for each bin the combinatorial background has been sub- 
tracted, as estimated from the corresponding toes dis- 
tribution, and the number of peaking background events 



has been estimated using fully simulated Monte Carlo 
events. The result of the fit is shown in Fig. [2 The frac- 
tion of B°^D-p'+(U50) and B°^D~ (ir+ir°) nr events 
in the mass window 620 < m W7r o < 920 MeV/c 2 is found 
to be smaller than 0.02 at 90% confidence level (C.L.) 
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FIG. 2: m w ^o distribution for the combinatorial-background 
subtracted B— >D ± 7r T 7r° sample, containing 16214 events. 
The solid line is the fit projection, consisting of the three inter- 
fering components described in the text and an t/ies peaking 
background contribution, indicated with the dashed line. 



The proper time interval At between the two B de- 
cays is calculated from the measured separation Az, be- 
tween the -B rec and B tag decay points along the beam 
direction. We determine the B rcc decay point from its 
charged tracks. The B tag decay point is obtained by fit- 
ting tracks that do not belong to B ICC , imposing con- 
straints from the B rcc momentum and the beam-spot lo- 
cation. We accept events with calculated At uncertainty 
of less than 2.5 ps and |At| < 20 ps. The average At 
resolution is approximately 1.1 ps. We use multivari- 
ate algorithms that identify signatures in the B tag decay 
products to determine ("tag") the flavor to be either a 
B° or a B° 0. Primary leptons from semileptonic B 
decays are selected from identified electrons and muons 
and from isolated energetic tracks. The charges of iden- 
tified kaons and soft pions from D* + decays are also used 
to extract flavor information. Each event with an esti- 
mated mistag probability less than 45% is assigned to 
one of six hierarchical, mutually exclusive tagging cate- 
gories. The lepton tagging category contains events with 
an identified lepton, while other events are divided into 
categories based on their estimated mistag probability. 
The effective efficiency of the tagging algorithm, defined 
as Q — Si ej(l — 2wi) 2 , where Ci and Wi are the efficiency 
and the mistag probability, respectively, for category i, 
is 30.1 ±0.5%. 

Since the expected CP asymmetry in the selected B 
decays is small, this measurement is sensitive to the in- 
terference between the b~^u and b— >c amplitudes in the 



decay of B t&g . To account for this "tagside interfer- 
ence", we use a parametrization different from Eq. 
which is described in Ref. 0] and summarized here. 
For each tagging category i, independent of the decay 
mode fi G {Dtt, D*n 7 Dp}, the tagside interference is 
parametrized in terms of the effective parameters r[ and 
<5j'. Neglecting terms of order (r^) 2 and (r^) 2 , the At 
distributions are written as 



ft'^^At) 



;-|At|/T 

4r 



[1 T K T - v4) 



sin(AmdAt) =p rj cos(AmdAt)] , (4) 



where, in the Standard Model, 

a* 1 = 2r^sm(2/3+j)cos5 fi , 

h = 2r-sin(2/3+7)cos<5-, 

cf = 2cos(2/3+ 7 )(r' i sm<5' J -r^sm^) 



(5) 



Semileptonic B decays do not have a doubly-CKM- 
suppressed amplitude contribution, and hence 



lep 



0. 



In the following, we quote results for the six a M and c(^ p 
parameters, which are independent of the unknown pa- 
rameters r- and 8[. The other hi and cf parameters de- 
pend on r[ and 5[, and do not contribute to the interpre- 
tation of the result in terms of sin(2/3+7). Note that all 
tagging categories contribute to the measurement of the 
a M parameters. 

An unbinned maximum-likelihood fit is applied to the 
At distribution of the selected B candidates in the AE 
signal region. The whole raES range is used to deter- 
mine the signal probability of each event on the basis of 
the Argus plus Gaussian fit described previously. The 
effect of finite At resolution is described by convolut- 
ing Eq. 0] with a resolution function composed of three 
Gaussian distributions. Incorrect tagging dilutes the pa- 
rameters a M , cf , and the coefficient of cos(AmdAt) by 
a factor D. L — 1 — 2wi Q. The parameters of the reso- 
lution function and those associated with flavor tagging 
are determined from the fit to the data and are consis- 
tent with previous BABAR analyses 0]. The At distri- 
bution of the combinatorial background is parametrized 
using two empirical components: a prompt component 
with zero lifetime and a component with an effective life- 
time. The components are convoluted with the sum of 
two Gaussians, and the resolution parameters of the two 
Gaussians, including the effective dilution parameters, 
the effective lifetime, and the relative fraction of the two 
components, are determined from the fit to the data. The 
peaking background coming from B ± mesons is mod- 
eled by an exponential with the B ± lifetime. Its relative 
fraction is fixed to the value estimated from simulations. 
The resolution function is the same as the signal resolu- 
tion, while the dilution parameters arc fixed to the values 
obtained from a B + control sample. The peaking back- 
grounds from B° mesons, whose amounts are also fixed to 
the value estimated using simulation, are modeled with a 



likelihood similar to the signal likelihood, but without CP 
violation (all the a, b, c parameters set to zero). Possible 
CP violation in this background is taken into account in 
the evaluation of the systematic uncertainties. The res- 
olution and the dilution parameters are the same as for 
the signal. 

From the unbinncd maximum likelihood fit we obtain: 



a D « = 

Dtt _ 

c lcp — 

a D ^ = 

ZTtt _ 

c lep — 

a D " = 

Dp _ 

c ic P ~ 



-0.010 ± 0.023 ±0.007 
-0.033 ±0.042 ±0.012 

-0.040 ±0.023 ±0.010 
0.049 ±0.042 ±0.015 
-0.024 ±0.031 ±0.009 
-0.098 ±0.055 ±0.018 



(6) 



where the first quoted error is statistical and the second is 
systematic. The largest correlation with any linear com- 
bination of other fit parameters is about 20% and 30% 
for the and the c(^ p parameters, respectively. Fig- 
ure |2| shows the fitted At distributions for events from 
the lepton tagging category, which has the lowest level of 
backgound and mistag probability. The various contri- 
butions to the systematic uncertainties of the a and ci ep 
parameters are shown in Table ITT1 



TABLE II: Systematic uncertainties on the a and ci ep param- 
eters (in units of 10~ 2 ). 



B { ' mode 


D ± 




D* ± -k t 


D* 


P T 


Source 




CTc 


£>~a 


o" c 




a c 


Vertexing (a At) 


0.37 


0.64 


0.80 


1.14 


0.47 


1.15 


Fit (cr flt ) 


0.51 


0.95 


0.52 


0.99 


0.75 


1.34 


Model (<r mod ) 


0.12 


0.13 


0.12 


0.13 


0.01 


0.18 


Tagging (utag) 


0.07 


0.16 


0.11 


0.14 


0.06 


0.12 


Background (obkg) 


0.13 


0.10 


0.10 


0.09 


0.28 


0.29 


m vn o Dependence (a\) 










0.16 


0.16 


Total (<7tot) 


0.66 


1.17 


0.97 


1.53 


0.94 


1.81 



The impact of a possible systematic mismeasurement 
of At (a At) has been estimated by comparing differ- 
ent parameterizations of the resolution function, varying 
the position of the beam spot and the absolute z scale 
within their uncertainties, and loosening and tightening 
the quality criteria on the reconstructed decay points. 
We also estimate the impact of the uncertainties on the 
alignment of the silicon vertex tracker (SVT) by repeat- 
ing the measurement using simulated events, with the 
SVT intentionally misaligned. For the systematic uncer- 
tainty of the fit (cat), we quote the upper limit on the bias 
on the a M and parameters, as estimated from samples 
of fully-simulated events. The model error (<T mo d) con- 
tains the uncertainty on the B° lifetime and Arrid, varied 
by the uncertainties on the world averages Q and also by 
allowing them to vary in the fit. The tagging error (otag) 
is estimated considering possible differences in tagging 
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FIG. 3: Distributions of At for the B°^D ± n T (a-d), 
B°^D* ± -k^ (e-h), and B°^D ± p =F (i-1) candidates tagged 
with leptons, split by B tagging flavor and reconstructed fi- 
nal state. The solid lines are fit projections. The background 
contributions are represented by the dashed curves. 



efficiency between B° and B°, different mistag fractions 
for the decay modes Dtt, D*tt, Dp, and different At res- 
olutions for correctly and incorrectly tagged events. We 
also account for uncertainties in the background (obkg) 
by varying the effective lifetimes, dilutions, wies shape 



parameters, signal fractions, and background CP asym- 
metry. A possible dependence of a Dp (c^) on the irir 
invariant mass (<J\) is estimated from the limit on the 
fraction of B°^D- p'+ (1450) and B°^D~ (ir + ir°) nr in 
the p mass window. 

As a cross-check, we perform the same fits on a sample 
of 6843 B~^D^°n~ candidates, where, as expected, 
we find no CP asymmetries. We combine our results 
with the result obtained on the partially reconstructed 
£?— >Z)* ± 7r f: sample [iH | and use a frequentist method de- 
scribed in Ref. jl(J to set a constraint on 2(3 + 7. The 
confidence level as a function of | sin(2/3+7)| is shown in 
FigureH We set the lower limits | sin(2/3+7)| >0.64 (0.40) 
at 68% (90%) C.L. 
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FIG. 4: Frequentist confidence level as a function of | sin(2/3+ 
7)|, obtained when combining our result with the result ob- 
tained on partially reconstructed B— >D* ± ir^ decays fiof . The 
horizontal lines show the 68% (top) and 90% C.L. (bottom). 

In conclusion, we have studied the time-dependent CP- 
violating asymmetries in fully reconstructed B a — ^D^tt^ 1 , 
B ^D* ± tt^, and B°^D ± pT decays in a sample of 232 
million T(4S)~^BB decays, and have measured the CP- 
violating parameters listed in Eq. We interpret the re- 
sult in terms of sin(2/3+j) and find that | sin(2/3+7)| >0.64 
(0.40) at 68% (90%) C.L. These results are consistent 
with and supersede our previous measurement. 
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